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1*.  Abstract 

Structural  Dynamic  Research  Cooperation  (SDRC)  obtained  FM 
magnetic  tape  from  Teledyne  Engineering  Services  containing  specified 
stress  records  from  three  Great  Lakes  vessels:  M/V  STEWART  J.  CORT, 
M/V  ROGER  BLOUGH,  and  S/S  EDWARD  L.  RYERSON.  Under  USCG  contract, 
SDRC  was  tasked  to  determine  the  damping  value  of  the  intervals  from 
these  records.  The  intervals  chosen  represented  periods  of  springing 
aboard  the  vessels. 

To  determine  a  relative  "feel"  for  the  accuracy  of  the  results 
obtained  by  SDRC,  LT  M.  D.  NOLL,  U.S.C.G.  Office  of  Research  and 
Development,  performed  an  evaluation  which  consisted 
of  a  manual  calculation  of  damping  and  logarithmic  decrement  from 
associated  quick-looks  and  a  comparison  of  these  with  the  results 
obtained  by  SDRC.  Although  the  comparison  does  not  show  complete 
agreement,  this  evaluation  does  conclude  that  the  SDRC  analysis  did 
provide  reasonable  values  and  may  be  valuable  for  future  data 
reduction  and  analysis. 
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Introduction 


During  the  past  two  decades,  vessel  designs  have  made  significant 
advances  which  include  increased  speeds,  reduced  hydrodynamic  resistance, 
increased  pay  loads,  and  increased  longitudinal  strength.  Oeat  Lake  ore 
carriers  have  also  made  significant  advances  but  with  a  special 
restriction  not  felt  by  ocean-going  vessels.  The  Oeat  Lake  ore  carrier 
has  a  max  inrun  beam  (B)  and  hull  depth  (D)  restriction  due  to  physical 
limitations  of  the  locks  at  Sault  Ste.  Marie,  Michigan.  However,  these 
Oeat  Lake  vessels  have  increased  in  length  to  1000  feet  with  an 
appropriate  lengthening  of  the  locks  by  the  U.S.  Amy  Corps  of  Bngineers. 

Of  interest  to  the  U.S.  Coast  Guard,  the  American  Bureau  of  Shipping 
(ABS)  and  the  shipping  industry  is  the  longitudinal  strength  required  of 
these  long,  shallow  ore  ships  exceeding  750'  in  length.  As  shown  in 
Table  1,  the  L/D  ratio  has  increased  from  18.7  to  22.23  giving  rise  to 
more  flexible  ships.  These  vessels  experience  a  phenomenon  known  as 
"springing",  a  2-node  vertical  oscillation  of  the  hull  structure,  caused 
by  the  wave  encounter  frequency  matching  and  being  in  phase  with  the 
natural  frequency  of  the  hull  girder.  The  addition  of  springing  to  the 
still  water  and  wave  induced  bending  mcments  may  cause  large 
cross-sectional  dynamic  stress  fluctuations.  These  increased  stresses 
and  associated  longitudinal  bending  mcments  must  be  accounted  for  in  the 
search  for  a  safe  and  effective  longitudinal  strength  standard. 

One  illusive  parameter  in  the  study  of  springing  has  been  the  damping 
factor  (  C  )•  System  damping  is  a  combination  of  the  hydrodynamic 
damping,  structural  damping,  cargo  damping  and  the  out-of-phase  seaway 
damping.  In  order  to  determine  the  damping  factor  for  Oeat  Lake  ore 
carriers  in  general,  stress  recordings  from  strain  gauges  mounted  on 
three  ore  carriers  have  been  analyzed  in  the  effort  to  determine  the 
damping  factor.  The  three  vessels  used  were  the  M/V  STEWART  J.  CCRT,  M/V 
ROGHl  BL0UCH,  and  S/S  EDWARD  L.  RYERSCN.  Table  1  lists  the  comparative 
characteristics  for  these  vessels.  The  stress  records  frcm  these  vessels 
have  been  recorded  and  analyzed  with  the  results  available  in  References 
1,2,  and  3. 
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TABLE  1  -  Vessel  Character istics 


CCRT 

BIOUCH 

RYERSGN 

LCA  (ft) 

998.0 

833 

730 

B  (ft) 

105.164 

105 

75 

L/B 

9.49 

7.9 

9.73 

D  (ft) 

44.9 

39.2 

39.0 

LA) 

22.23 

21.3 

18.7 

BA) 

2.34 

2.68 

1.92 

G.T.  (tons) 

32930 

22041 

12170 

N.T.  (tons) 

29918 

14114 

7637 

H.P. 

14000 

14000 

9000 

Yr.  Built 

1971 

1971 

1960 

Owner 

Beth.  Steel 

U.S.  Steel 

Inland  Steel 

Background 

In  June  1977,  Teledyne  Engineering  Services  (formerly  Teledyne 
Materials  Research)  of  Waltham,  Massachusetts  was  contracted  by  the  U.  S. 
Coast  Qiard  to  transfer  specified  intervals  of  analog  tape  stress  records 
frcm  these  three  vessels  onto  FM  magnetic  tape.  These  intervals  were 
chosen  to  represent  high,  medium  and  low  springing  stresses  experienced 
by  the  vessels  during  both  the  loaded  and  ballasted  condition.  The 
intervals  were  chosen  by  the  Coast  GLiard  frcm  the  computer  tables  listed 
in  Reference  1,2,  and  3.  Tables  2-4  (Colums  1,  2,  3  and  14)  provide  the 
listing  of  the  intervals  and  maxinum  peak-to-through  (P-T)  stresses  which 
Teledyne  was  contracted  to  copy  over  onto  magnetic  tape.  Tape  Ninber  and 
interval  number  correspond  to  computer  tables  in  these  references.  (For 
example:  Table  2,  Tape  No.  ”18",  Interval  No.  "20"  corresponds  to  same 
Tape  Number  and  Interval  Number  in  Reference  1,  Appendix  B,  page  48). 

Upon  receipt  of  theFM  magnetic  tapes  and  the  quick-looks  frcm 
Teledyne,  Structural  Dynamics  Research  Corporation  (SDRC)  of  Cincinnati, 
Ohio  was  contracted  to  determine  the  range  of  damping  values  for  these 
specific  intervals.  SDRC  took  possession  of  the  tapes  and  the 
quick-looks  for  their  damping  analysis. 

8EFC  Determination 


For  each  specific  interval,  SCRC  supplied  a  table  of  estimated  roots 
and  an  associated  damping  power  spectra.  The  tables  of  modal  parameters 
and  the  data  plots  are  connected  by  a  FRBCJIESP-BCDE  number  located  at  the 
top  of  the  table  and  at  the  bottom  of  the  plot  (i.e.  for  pages  A-l  and 
A-2,  this  number  is  IX  +  2Y+).  The  BCDE  nurber  before  the  V  refers  to 
the  "ID"  colum  in  the  table  of  calculated  damping  values  on  pages  A-3 
and  A- 4.  This  allows  the  correct  correlation  with  appropriate  table  and 
plot.  Appendix  F  provides  a  complete  listing  of  all  damping  power 
spectres  submitted  by  SCRC. 
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Figure  6  (a,b,e)  shows  a  plot  of  amplitude  versus  frequency  for  the 
CXHT,  BLOUCH,  and  RYHISQN  respectively,  which  illustrates  how  the  max 
amplitude  under  this  procedure  can  be  used  to  choose  the  appropriate 
damping  factor.  Notice  that  wtien  all  data  is  plotted  as  amplitude  versus 
frequency,  two  peaks  appear.  This  may  be  misleading  since  the  peak  of 
interest  occurs  at  the  known  natural  frequency  for  each  vessel. 

As  mentioned  before,  the  damping  values  calculated  by  SERC  represents 
the  total  system  damping  which  includes  hydrodynamic,  structural  and 
cargo  damping.  The  representative  system  damping  factor  and  frequency 
for  the  specific  interval  was  chosen  by  picking  the  damping  factor  and 
frequency  that  corresponded  to  the  max  amplitude  of  the  associated  power 
spectrun.  For  example:  the  damping  factor  of  0.01247  chosen  from  page 
A-l  corresponds  to  a  frequency  of  .3455  Hz.  which  corresponds  to  the  max 
amplitude  of  0.4804E-03  shown  in  the  table  of  estimated  roots  (pg.  A-3) 
and  plotted  on  page  A-2  (As  mentioned  in  the  cover  letter  from  SCRC 
(Appendix  B),  the  amplitude  is  in  error  by  a  constant  factor). 

In  order  to  portray  the  complete  picture  at  a  glance,  Tables  2^  3  and 
4  were  expanded  to  include  other  pertinent  data  such  as  ship  course  and 
speed,  wind  direction  and  speed,  wave  direction  and  height,  max  P-T 
stress,  and  etc.  This  data,  extracted  from  References  1,2,  and  3,  is 
listed  in  Table  2,  3,  and  4  to  correspond  to  the  OCRT,  BIOUGH,  and 
KTERSON,  respectively.  Parentheses  used  in  the  Tables  indicate  that  the 
data  was  extracted  fran  the  last  recorded  entry  in  that  colum  from  the 
original  tables  in  References  1,  2,  and  3. 
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Method  of  Analysis 

The  system  of  analysis  utilized  by  SERC  is  known  as  MODEMS  (MQDal 
Analysis  and  Modeling  System).  The  technique  used  to  determine  the 
desired  nodal  parameter  estimations  was  the  mil ti -degree -of -freedom 
(MDOF)  curve  fit  to  the  frequency  response  data  over  a  specific  frequency 
range,  fsee  Appendix  E  for  a  brief  description  of  a  single  d.o.f.  system 
used  by  ABS) 

The  algorithm  used  by  SCRC  is  explained  in  Appendix  B  of  this 
report.  This  description  gives  the  theoretical  background  for  the 
curvefit  algorithm  used  to  extract  modal  parameters  from  the  response 
spectra.  This  algorithm  will  work  on  spectra  with  no  phase  information 
but  results  in  an  analytical  function  of  frequency  whose  calculated  phase 
angle  is  essentially  meaningless.  Even  though  mode  shape  coefficients 
would  also  be  meaningless,  the  mode  shape  evident  in  each  graph  is 
damping  dependent  and  therefore  allows  the  determination  of  the  damping 
factor  despite  the  lack  of  the  phase  information. 

The  basic  concept  behind  the  extraction  of  the  damping  values  centers 
on  the  assumption  that  the  operating  data  consists  of  a  system's  response 
to  a  fairly  randan  forcing  function  (fairly  flat  frequency  of  input). 

The  measured  response  can  be  strain,  displacement,  velocity,  or 
acceleration.  The  first  two  are  used  directly  and  the  latter  two  are 
converted  to  displacement  data  first.  Since  strain  relates  to 
displaconent  directly,  the  MDOF  technique  calculates  Af|<  (page  7, 
Appendix  B)  fran  a  displacement/force  function.  Since  SERC  assured  the 
force  applied  was  a  constant  function  of  frequency,  using  F(w)  =1.0  lbf 
for  division  in  X/F  did  not  affect  the  overall  shape  of  the  response 
curves. 


Using  this  method  allowed  SERC  to  establish  the  tables  and 
associated  graphs  for  each  interval  specified  by  the  U.  S.  Coast  Guard. 
The  damping  factor  and  frequency  range  chosen  for  each  interval 
corresponds  to  the  max  amplitude  for  the  frequency  range  investigated, 
specifically  0.01  to  1.0  Hz. 

Analysis  of  SERC  Results 

The  range  of  values  supplied  by  SERC  still  raises  the  question  of  the 
accuracy  of  their  results.  Determining  the  accuracy  may  be  difficult  but 
determining  whether  or  not  they  are  in  the  "ballpark"  would  be  highly 
beneficial . 


The  damping  value  (  C  )  calculated  by  SERC  is  the  fraction  of  critical 
damping  where:  C  =  c 

®c 

and:  c  =  damping  coefficient,  tons-sec/in 

cc  =  critical  damping  coefficient,  Ib-sec/in 
=  arwn  =  2  "vnk 

wn  =  natural  frequency  =  Vk/m 
k  =  stiffness  factor,  tons/in 
m  =  mass,  tons-sec^/in 


Then  C  represents  a- certain  ratio  of  e  to  cc  to  arrive  at  the 
damping  coefficient  c  where  cc  *  C  =  c. 

In  order  to  determine  a  "feel”  for  the  approximate  magnitude  of  the 
damping  values  for  the  three  vessels,  a  mechanical  method  was  used  to 
determine  the  different  springing  parameters.  Tins  method  utilizes  the 
filtered  records  fron  existing  stress  records.  It  allows  the 
determination  of  the  approximate  damping  values  from  mechanical 
measurements  taken  directly  from  stress  records. 

Suppose  that  a  free  transient  has  the  following  form: 


Wiere: 


Figure  1  -  Example  of  a  Free  Transient 

X(t)  =  Amplitude  at  time  t,  psi 
Xq  =  Initial  amplitude,  psi,  at  t=0 
Xj  =  Amplitude  after  initial  oscillation,  psi 
X,,  =  Amplitude  after  n  oscillations 
n  =  Nlrrber  of  oscillations 
tn  =  Time  for  n  oscillations  to  occur,  sec 
t(j  =  Damped  period,  sec 
w<j  =  Damped  Frequency 
wn  =  Natural  Frequency 
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The  free  transient  response  for  this  syston  shown  in  Figure  1  is: 


Xo(t)  =  Xq  *  e~Cte>|t  *(  sin (  Vl  -g2  Wnt  +  +  )) 

v  i  - 

With  the  initial  conditions  being: 


X(0)  =  Xo 
X(0)  =  0 


These  initial  conditions  are  used  in  the  mechanical  analysis  because 
the  actual  forces  causing  a  peak  condition  is  of  no  concern  to  this 
analysis.  The  area  of  interest  is  the  decay  period  of  the  record.  Chce 
the  peak  springing  stress  is  reached  on  a  specific  interval,  the  analysis 
is  measuring  the  rate  of  decay  as  the  oscillations  decrease  in  anplitude. 


Fran  the  springing  record,  Xq,  X„,  and  n  may  be  measured.  Fran 
these  values  the  damping  can  be  determined  as  follows: 


XO 


Xq(0) 


Xo(0) 


Xn  X„(tn) 

and  rearranging: 


Xo(  n  *  td) 
1 


e  *  2irn 


•  In 


XQ 


Yl  -  cz  2 it  n 

When  c  «  1.0,  then  Vl  -j2  =  1.0.  Therefore: 

1  Xo 


Xn 


t  = 


*  In 


(1) 


2  irn 


Xn 


Equation  1  is  the  formula  used  to  calculate  the  danping  values  found 
in  Table  5.  Appendix  C  contains  the  stress  records  from  which  the 
appropriate  measurements  were  made  in  order  to  calculate  damping  via  the 
mechanical  method.  From  the  same  records  the  damped  frequency  (  wd), 
the  natural  frequency  (  wv,),  the  damped  period  (td),  and  the 
logarithnic  deeronent  (fi  )  may  be  calculated  via  the  following  formulas: 

wh  =  n  Hz.  =  2  ffn  radians  (2) 

tn  tn  sec. 
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1 


* 


Hz. 


(3) 


r 


wd 


Note:  as  per  conditions  for  Equation  1,  w<j  *  w^, 

*n 

t<j  = _  sec. 

n 


(4) 


The  term  In  Xq/X,,  from  Equation  1  is  often  referred  to  as  the 
logarithmic  decrement  and  is  defined  as  the  natural  logarithm  of  any  two 
successive  amplitudes.  Since  Table  5  lists  the  peaks  amplitudes,  Xq, 
and  the  amplitude  after  n  oscillations,  Xf,,  the  equation  for  the 
logarithmic  decrement  is  modified  to  be: 

6  =  _L  •  1"  Jfo.  (5) 

n  Xn 

These  parameters  are  calculated  from  the  measurements  taken  from 
Appendix  C  and  the  results  are  incorporated  into  Table  5.  The  areas  of 
measurement  are  enclosed  by  double  rectangles  and  are  nimbered  to 
correspond  to  the  miner ical  order  in  Table  5. 

Logarithmic  decrement  may  also  be  determined  if  the  damping  value  is 
already  known.  The  equation  for  4  can  be  written  in  the  simplified  form: 

6  =  2*C  =  2  it  c  (6) 

Vl  -  C2 

Equation  6  is  an  alternate  method  for  determining  fl  when  C  is  known  and 
the  stress  records  are  not  available  for  mechanical  measurements  as  performed 
for  Table  5. 

Comparison  of  Results 

In  order  to  provide  a  means  of  comparison  between  the  SDRC  results  and  the 
mechanical  results,  a  program  was  written  in  BASIC  to  calculate  the  <5  ,  t<j, 
the  bending  moment  (£M)  from  the  data  presented  in  Tables  2,3,  and  4.  A  flow 
chart,  the  program  data  and  the  results  are  shown  in  Appendix  D. 


i 


r 


Table  S  -  Damping  Parameters  Calculated  from  Mechanical  Measurements 


■■ 

ftl 

H 

■1 

H 

*0 

psi 

*n 

psi 

n 

ey 

(sec) 

C 

Wd 

Hz. 

wn 

Hz. 

6 

Ref.  • 

1 

8206 

2500 

7 

20.2 

.0207 

.35 

.35 

2.89 

.129 

C1:R9:14 

2327 

517 

5 

13.4 

.048 

.37 

.37 

2.68 

.30 

Cl  " 

H 

4913 

2586 

6 

17 

.017 

.35 

.35 

2.83 

.107 

C2  " 

6465 

1552 

5 

14.6 

.0434 

.34 

.34 

2.9 

.285 

C2  " 

5 

7758 

1034 

8 

24 

.040 

.33 

.33 

3.0 

.251 

C3  " 

6 

7500 

3879 

4 

12.2 

.026 

.33 

.33 

3.05 

.163 

C3  " 

7 

7241 

1243 

9 

26.6 

.0311 

.34 

.34 

2.96 

.191 

C3  " 

8 

6350 

2533 

5 

17.6 

.0293 

.28 

.28 

3.52 

.184 

C4  " 

9 

11200 

1867 

18 

52.7 

.0158 

.34 

.34 

2.93 

TO 

C5:Q-L 

10 

1917 

800 

6 

17.6 

.0232 

.34 

.34 

2.93 

.146 

C6:Q-L 

BIOUCH 

11 

6600 

1000 

12 

30.0 

.025 

B 

B 

2.5 

.025 

C7:R2 

12 

9250 

2434 

7 

15 

.0304 

H 

2.14 

.191 

C7  " 

13 

5780 

2720 

6 

16 

.02 

2.67 

.126 

C7  " 

RYERSCN 

14 

(.28) 

7 

12.5 

.0289 

.56 

.56 

1.79 

C8:R3 

15 

(.47) 

— 

5 

9.0 

.0240 

.56 

.56 

1.80 

y 

08  " 

•Reference  Code  C  -  Appendix  C  of  this  report 

Cl  -  pg.  1  of  Appendix  C  of  this  report 
R9  -  Reference  9  (page  19) 

14  -  pg.  14  of  Reference 
Q-L  -  Quick-look 

••NUrbers  in  parenthesis  is  %  with  Xq  being 
unit  1  and  Xf,  being  a  fraction  of  Xq. 
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The  locator  code  in  the  ccmputer  output  refers  to  Tables  2,3,  and  4. 
For  exanple:  the  first  locator  code  in  the  computer  output  (page  D-4)  is 
"CHSL*18*20".  This  refers  to  Table  2  meaning  "CCRT-High  Springing  Load, 
Tape  #18,  Interval  #20".  In  this  way  a  direct  connection  nray  be  made 
between  Tables  2,3,  and  4  and  the  ccmputer  output.  The  equations  used 
for  the  calculations  of  the  variables  in  the  program  are  located  on  page 
D-2. 


Using  the  output  from  the  computer,  Tables  1  -  5,  and  the  SBC 
results,  a  comparison  was  performed  to  see  how  close  all  inputs  agree. 
Figure  2  is  a  plot  of  6  versus  C  .  TTie  "Key"  differentiates  between 
the  SERC  results  and  the  mechanical  results.  Since  t  is  small  (i.e.<< 
1.0),  it  was  expected  that  the  data  points  would  lie  on  a  straight  line 
since  s  is  directly  proportional  to  t  as  per  Equation  6.  The  nunbered 
points  on  Figure  2  correspond  to  the  sequential  ntmbering  from  Table  5. 
Gfenerally  the  mechanical  results  plotted  linearly  with  the  SRC  results 
but  gave  higher  values  of  C  .  This  could  possibly  be  explained  by  the 
lack  of  more  information  on  the  actual  sea  state  existing  at  the  time  of 
stress  measurements.  Table  6  below  gives  the  relative  comparisons 
between  the  two  results. 

Table  6  -  Comparison  of  SRC  Results  with  Mechanical  Results 


SRC 

Mechanical 

Range  of  5 

Avg. 

Range  of  C 

Avg. 

SRC  avg  *  100 
Mech.  avg. 

CERT 

.00611-. 0246 

.014 

.017-. 048 

.029 

48  % 

BLOUCH 

.0103-. 0248 

.01916 

.020-. 0304 

.0251 

76  % 

RYFRSON 

.00337-. 00745 

.00480 

.024-. 0289 

.0265 

18  % 

As  can  be  seen  from  Table  6,  the  differences  in  the  calculated  damping 
factors  are  quite  pronounced,  especially  for  the  OCRT  and  the  RYER90N.  At 
first  inspection,  this  would  tend  to  invalidate  the  SRC  results.  However, 
despite  the  observed  differences  in  Table  6  and  Figure  2,  an  examination  of 
frequency  versus  damping  yields  an  interesting  results. 

Figure  3  is  plot  of  w  versus  c  using  the  results  from  SCRC.  The  data 
points  for  each  vessel  are  enclosed  within  a  dashed  line.  From  Table  1  it  may 
be  seen  that  the  CERT  has  a  higher  L/D  ratio  (22.2)  than  the  other  two  vessels 
and  this  would  tend  to  make  the  CERT  more  "flexible"  than  the  BLOUCH  or  the 
RYFR80N.  This  in  turn  would  cause  the  CERT  to  have  a  lower  springing 
frequency  than  the  BLOUCH  or  the  RYERSON.  Since  the  RYHtSON  has  the  lower 
LA)  ratio  of  18.7,  it  would  tend  to  be  a  stiffer  ship  and  show  a  higher 
springing  frequency.  These  observations  are  shewn  in  Figure  3.  The 
approximate  springing  frequency  for  the  CERT  is  0.32  Hz.  despite  the  loading 
condition  and  despite  the  stress  levels  seen.  By  reviewing  the  original 
stress  levels,  the  Low  Springing  Ballast  case  for  the  CURT  is  approximately 
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one-third  the  levels  in  the  High  Springing  Load  ease.  (See  Table  2) 

This  is  an  important  observation  illustrating  that  exciters  that  would  be 
used  for  establishing  the  fundamental  mode  of  hull  vibration  (springing) 
would  not  have  to  generate  peak  stress  levels.  The  box  girder  of  the 
vessel  would  still  spring  at  approximately  0.32  Hz.  The  damping  would  be 
measured  more  accurately  since  the  out-of-phase  seaway  damping  would 
essentially  be  eliminated. 

Figure  4  is  a  plot  of  damping  versus  frequency  as  calculated  in 
Table  5.  The  nurtoered  points  correspond  to  the  sequential  numbering  in 
Table  5.  The  enclosed  dashed  lines  represent  the  SERC  results  in  Figure 
3.  As  shown,  several  points  are  outside  the  representative  areas  but  the 
corresponding  points  do  correspond  to  the  same  frequency  as  depicted  in 
Figure  3.  Comparison  between  Figure  3  and  Figure  4  illustrates  this 
notable  point:  using  only  analog  stress  records  (digitized  by  Teledyne) 
during  periods  of  vessel  springing,  SERC  developed  a  list  of  damping 
factors  and  associated  frequencies  at  spr inging  which  matched  "fairly" 
well  with  a  mechanical  analysis  from  apparently  disassociated  springing 
records  (Appendix  C)  and  agrees  with  the  known  natural  frequency  of  each 
vessel . 

The  unknowns  in  both  analyses  are  (1)  the  components  of  the  total 
damping,  (2)  the  method  and  distribution  of  seaway  loading  on  the  vessel, 
(3)  the  out-of-phase/ in-phase  characteristics  of  the  seaway  with  the  hull 
response,  (4)  understanding  of  the  full  complexity  of  the  hull  structure, 
and  (5)  more  accurate  and  complete  description  of  the  sea  conditions 
(i.e.  wave  height,  wave  length,  wave  direction,  wind  parameters,  and 
actual  loading  conditions).  With  this  as  part  of  the  total 
consideration,  the  results  from  SEFC  are  very  encouraging  despite  the 
differences  calculated  in  Table  6. 

As  one  final  attempt  for  correlation  between  the  two  methods,  the 
program  in  Appendix  D  also  calculated  the  bending  mcment,  BA,  using  both 
the  P-T  stresses  of  each  vessel  from  Appendix  A,  but  also  using  the  max 
stress,  Xq,  fron  Table  5.  TTie  calculated  EM  was  plotted  versus  the 
damping  factor.  The  plot  is  shewn  as  Figure  5. 

TTie  values  of  BA  calculated  by  Appendix  D  are  shown  in  Column  6  of 
page  E>-4  of  Appendix  D.  The  values  for  BA  from  Table  5  are  shown  in 
Table  7.  The  EM  was  calculated  by: 


04  =  3A  »  o _  FT-Tons 

12  *  2000 


where: 

BA  =  total  bending  moment. 
5M  =  section  modulus,  in.3 
o  =  stress,  psi 
12  =  12  in. /ft. 

2000  =  2000  lb. /ton 


ft -tons 
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Table  7  -  Calculated  Bending  M^nents  from  Initial  Stress  of  Table  5 


Xo 

m 

1. 

6206 

PSI 

23833 

FT-TONS 

2. 

2327 

n 

8986 

II 

3. 

4913 

it 

18868 

If 

4. 

6465 

it 

24828 

II 

5. 

7758 

11 

29794 

If 

6. 

7500 

r? 

28800 

1! 

7. 

7241 

if 

27808 

If 

8. 

6350 

if 

24386 

II 

9. 

11200 

if 

43012 

II 

10. 

1917 

n 

7361 

II 

11. 

6600 

r» 

16500 

11 

12. 

9250 

it 

23125 

If 

13. 

5780 

r? 

14450 

II 

The  SM  used  for  each  vessel  is  shown  on  page  D-2.  Since  the  HYStSON 
worked  by  relative  values  of  Xq  and  Xn,  no  IM  was  calculated.  TTie 
calculated  values  for  each  vessel  as  per  the  SERC  results  were  boxed  in. 
The  B4’ s  calculated  in  Table  7  were  overlayed  with  the  numbered  points 
corresponding  to  the  sequential  nurfcering  in  Table  5  and  in  Table  7. 

TTere  appears  to  be  "general"  agreement  but  there  appears  to  be  no 
conclusions  that  can  be  derived  frcm  this  plot.  The  plotted  points  would 
be  in  better  agreement  if  the  damping  factors  were  in  better  agreement. 
Therefore,  the  points  for  the  BIOUCH  show  good  agreement  since  the  SERC 
results  are  76%  of  the  mechanical  results,  (see  Table  6)  Che  final 
carment  on  Table  6:  just  because  the  results  show  better  correlation  of 
the  SERC  results  with  the  mechanical  results  for  the  M/V  BLOUCH  does  not 
mean  that  a  better  analysis  was  performed  on  the  BIOLJCH  or  that  better 
records  were  achieved  or  that  springing  was  better  identified.  The  only 
conclusion  that  can  be  drawn  as  result  of  Table  6  is  that  the  comparative 
springing  records  in  Appendix  C  may  have  more  closely  resembled  the 
BIOUCH  records  used  by  SERC.  The  sea  conditions  may  have  been  closer  to 
being  the  same.  There  is  nothing  to  say  that  if  different  springing 
records  had  been  chosen  for  conparison,  the  GCRT  or  RYERSON  results  may 
have  proven  to  give  the  best  comparison  with  the  SERC  results. 

Conclusion 


Considering  the  information  that  SHC  had  and  did  not  have,  the 
results  of  their  analysis  is  very  encouraging.  The  damping  values 
provided  are  a  combination  of  structural  damping,  hydrodynamic  damping, 
cargo  damping,  and  out -of  phase  seaway  damping.  The  data  provided  by 
SERC  may  be  used  for  further  analysis.  In  order  to  arrive  at  the  true 
structural  damping  within  a  vessel,  an  exciter  capable  of  producing 
springingmay  be  needed  in  order  to  eliminate  so  many  of  the  unknowns 
which  have  complicated  this  analysis.  The  sign  and  operation  of  such  a 
device  presents  its  own  difficulties  in  the  form  of  unknowns, 
assumptions,  and  feasibility. 
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Figure  5  -  Bending  Moment  vs.  Damping  Factor 
Comparison  of  SDRC  and  Mechanical 
Results 
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Figurt  6(b)  -  Frequency  vs.  Amplitude  -  BLOUGH 
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Figure  6(c)  -  Frequency  vs.  Amplitude  -  RYERSON 
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APPENDIX  A 

Sample  Results  from  SDRC 
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APPENDIX  B 

^'^QPStical  Background  for  Curve  Fit 
Algorithm 


Structural  t'/rri'.iLs  Reseat ch  Corporation 
l»/29  D*  agon  Way 
Cincinnati,  Ohio  45227 
5  id- 2  7 2-1 100  TWX  810-461-2615 


August  15,  1977 


Lt.  M.  Noll 

Commandant  (G-DSA/TP44) 

U.S.  Coast  Guard 

Washington,  D.C. 

Subject:  Results  of  Springing  Damping  Analysis,  SDRC  Proposal  No.  6704-01 ;  Your  RFQ  No 
G-FCP-13/C22/77. 

Dear  Mark: 

Attached  you  will  find: 

1.  Table  of  calculated  damping  values. 

2.  Tabulation  of  modal  parameters  from  each  curve-fit 

3.  Plots  with  superposition  of  data  and  curve  fit 

4.  Theoretical  background  for  curve-fit  algorithm 

There  are  a  few  comments  I  would  like  to  make  concerning  the  interpretation  of  results: 

1.  Fhe  curve-fit  algorithm  used  on  this  data  to  extract  modal  parameters  is  briefly  described 
in  the  attachment  (4).  The  form  of  the  data  we  curve:fit  was  PSD,  of  course,  not  transfer 
functions.  However,  the  data  was  assumed  to  be  a  transfer  function  where  the  force  applied 
was  a  constant  function  of  frequency,  F(w)  =  1.0  Ibf  so  that  division  by  this  to  get  X/F 
(actually  strain/force,  which  is  not  an  important  difference  here)  will  not  affect  the  shape 
of  the  curves. 

2.  The  algorithm  we  use  to  curve  fit  transfer  functions  will  work  on  spectra  with  no  phase 
information.  The  result  is  simply  an  analytical  expression  as  a  function  of  frequency  whose 
amplitude  is  in  error  by  a  constant  factor,  and  of  course  the  calculated  phase  angle  is  meaning¬ 
less.  Thus  mode  shape  coefficients  calculated  ftom  a  power  spectrum  would  be  meaningless. 
The  damping  on  the  other  hand  does  not  suffer  from  the  lack  of  phase  angle.  As  the  super¬ 
imposed  ptots  show  the  overall  shape  of  the  curves,  which  ace  damping  dependent,  are  correct, 
but  the  plot  must  be  shifted  up  or  down  the  vertical  axis  to  obtain  the  correci  amplitude.  The 
constant  by  which  the  curve  is  in  error  depends  upon  the  frequency  of  the  most  dominant 
peaks,  so  not  all  the  plots  show  the  same  amplitude  error. 

3.  The  Hilbert  transform  is  used  to  calcualte  a  phase  angle  from  a  power  spectrum  as  a  function 
of  frequency.  The  inclusion  of  the  phase  angle  da'a  with  the  amplitude  data  in  the  curve  fit 
allows  exact  amplitude  correlation  between  transfer  functions  calcualted  from  the  curve  fit 
parameters  and  the  original  data  (see  Figures  1  and  2).  However,  as  the  enclosed  plots  show, 
the  visual  correlation  can  be  done  very  well  by  comparing  the  shapes  of  the  two  curves. 
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4.  The  reason  the  Hilbert  transform  was  not  utilized  in  this  project  is  that  hardware  problems 
with  the  computer  were  inhibiting  a  function  required  in  the  phase  angle  computation.  For 
expediency,  I  elected  to  proceed  with  the  data  reduction  rather  than  await  the  correction 
problem.  The  algorithm  for  curve-fitting  was  checked  out  on  other  data  prior  to  reduction 
of  the  Coast  Guard  data  to  verify  the  hardware  problem  was  not  interfering  with  its  opera¬ 
tion. 

If  you  need  more  information,  or  have  any  questions,  please  contact  me. 


Regards, 


if!..  M-'VJjSiU  ithjfa — 


R.  Gene  Smiley 
Project  Manager 
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APPENDIX  A  -  THEORETICAL  BACKGROUND 

The  theory  behind  nodal  analysis  via  frequency  response  functions 
can  be  exaoined  by  referring  to  the  equations  of  notion  for  an 
N  degree  of  freedom  system  with  viscous  danpingj 

0:][q]  +  [C][q]  +  [K][q]  =  [f]  (1) 

where 

[M]  =  mass  matrix 

[C]  =  viscous  damping  matrix 

£k3  =  stiffness  matrix 

£q]  =  time  history  of  the  displacement  of  system 
[f]  =  time  history  of  excitation  to  system 
=  time  history  of  velocity  of  system 
[q]  =  time  history  of  acceleration  of  system 
This  equation  is  inconvenient  to  handle  with  standard  methods  of 
eigenvalue  analysis  if  [c]  is  not  proportional  to  [.*.:]  or  [K] . 
However,  a  method  has  been  proposed  by  Dur.can  (1)*  which  reduces 
these  equations  to  a  standard  eigenvalue  form.  In  this  method 
combine  the  identity* 

DOCq]  -  DO  [q]  =  [0] 

with  Equation  1  to  obtain* 


f[0]  CM] 

ra]  + 

r-Daw 

[q] 

-  -t 

[0] 

l pa  [c] 

L[o][k] 

Cq] 

•  '  m 

Represent  this  equation  in  the  following  manner* 

M[y)  ♦  [3][y]  -  CO 

l 

•lumbers  in  parentheses  designate  references  at  end  of  s=c~ior 
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(3) 


In  order  to  find  the  solution  to  2o.ua t ion  3  for  the  case  of  harmonic 
inputs t  first  consider  the  solution  to  the  homogeneous  equation 
found  toy  letting  [z]  =  [0] 

[A][y]  +  [B][y]  -  [o]  (4) 

Seek  a  solution  of  the  form  [y]*=  £Y]est 
therefore  =  s  [Y]est 
Hence,  Equation  4,  becomes, 

sCA][y]  +  C3][y]  =  C°]  °r 

[Cb]  +  »[A]]  [y]  =  [0] 

This  set  of  equations  only  have  a  solution  if  the  determinant  of 
the  coefficient  matrix  is  zero, 
det  £[3]  +  s[A]J  =  C 

This  leads  to  a  set  of  2N  roots  or  eigenvalues  sq,  S2  •  •  •  •  S2n» 
which  satisfy  the  above  equation.  For  a  resonant  system,  these  eigen¬ 
values  will  occur  in  conjugate  pairs  (2).  Corresponding  to  each 
eigenvalue  sr,  there  exists  an  eigenvector  [v}'r]  having  21'i  components 
satisfying  the  following  equation! 

[C3]  +  SrCA]J  C-PT-Cc]  (5) 

In  the  case  where  the  eigenvalues  of  a  system  are  complex,  in.  which 
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case  they  occur  in  conjugate  pairs,  the  eigenvectors  will  be 
complex  and  will  also  occur  in  conjugate  pairs, 


The  above  eigenvectors  have  important  orthogonality  conditions 
which  can  be  easily  shown.  Consider  the  r*11  and  p*h  eigenvectors 
£  C  both  of  which  satisfy  Equation  5*  First  write 

Equation  5  for  the  r^*1  mode  and  premultiply  by  the  transposed 
vector  obtain i 

['fp]TCB][vpr]-3rC'J'p]TCA]['i:r]=r03  <«> 

Using  the  reversal  law  for  transposed  matrix  products  and  re¬ 
calling  that  [a]  and  £3]  are  symmetric  matrices,  transpose 
Equation  6  to  obtain: 

C?r]T[B]['fp]+3r[9r]T[A][f  P]«[o]  (?) 

Next  write  Equation  5  for  the  pth  mode  and  premultiply  by  £{,r]‘ 


[i,r]TCB][f  P]+Sp[  i  r]T[A][£  P]=L°]  (8)  ■ 

If  Equation  8  is  subtracted  from  Equation  ?,  the  result  isi 

<vV  C?r]TWCfp]  =  o 

If  eigenvalues  sr  and  sp  are  different,  the  following  orthogonal! 
property  relates  the  two  eigenvectors. 

CTr],WC4p]  -  o  (9) 

It  follows  that  these  vectors  are  also  orthogonal  with  respect 
to  matrix  £a] 

C}r]T[B]['?B]  =  0  (10) 
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Equations  9  and  10  are  important  orthogonality  conditions  which 
shows  that  the  2N  vectors  [SE'l  form  a  linearly  independent  set, 
•and  therefore  any  vector  in  2K  space  can  be  expressed  as  a 
linear  combination  of  these  2H  vectors.  Since  we  are  interested 
in  frequency  response  information,  let 

ouiHz]  e 3wt  an 

and  seek  a  solution  in  the  forms 

[y(t)J  =[Y]ej:ot  (i2) 

Substitute  Equation  11  and  Equation  12  into  Equation  3  and 
divide  by  eJ  to  obtain: 

jco[A][Y]  +  [B][Y]  =  [Z]  (13) 

Since  the  eigenvectors  defined  by  Equation  5  form  a  linearly 
independent  set  over  2N  space,  write  the  solution  to  Equation 
(13)  as  a  linear  combination  of  these  2N  vectors 


[Y]  * 


xr£Trl 


Substitute  Equation  14  into  Equation  13  and  multiply  by  [  'p?]'1’ 
to  obtain: 

2M  2Nr  - 


p]Tu]  vrro'r>»p]T[B 


r=l 


From  the  orthogonality  conditions,  we  obtain: 


i“V  V  •  V 


where 


ap  *  p]tCa]C-Pp] 

*p  *  [¥?]TL3][^p] 
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Hence,  we  can  solve  Equation  15  for  to  obtaim 

DggftXl 

**  j<Oap+  bp  (16) 

Substituting  Equation  16  into  Equation  14,  we  obtain: 

2N_  [^r]T[z]['fr] 

(17) 


m  = 


^=T  3^ar  + 

However,  from  Equation  7i  we  obtain 


br  +  sr  ar  =  ° 

_  “by 


Therefore,  Equation  (17)  can  be  '.written 
2N  C^r]T[Z][vpr] 

ra  -  ^  — —■  /  as) 

ar(jcj  -  sr) 

Frequently  the  complex  eigenvalues  sr  are  written  in  the 
following  form: 

*r"-'£rUr  *1  "r/1*--5,' 


where 

3. 


damning  ratio 


cor  =  undamped  natural  frequency 
In  terms  of  [Q]  and  [?] ,  Equation  18  becomes: 


JH-  [^r]T[g]  [fr] 


-  J«o|  Lj 

.  Q  J  ar 


(jw-  ^rdjr+j  ox^A-  ) 


Therefore  the  frequency  response  function  recorded  from  excitation 
applied  at  location  k  and  response  monitored  at  location  i  is: 
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(A-6) 


Hik= 


W  ?iT _ 

;z^‘ar(  j  CO  +  ^r±  j  Or/l-  3^) 


Since  the  eigenvalues  occur  in  conjugate  pairs,  Equation  19 
can  be  written  as 


Hik= 


jr  ?kr  ^ir  + 

r=l  ar(jco  +  Sr«Jr+j  ^ry 

¥/*  _ 

a?u^+  or  cjr  -j<jri/rr^) 


Equation  19  is  an  extremely  valuable  relationship  be  tv/ sen 
Frequency  Response  Functions  and  modal  characteristics.  It 
relates  motion  at  any  point  i  due  to  an  force  at  point  k. 
Notice  that  equation  19  implies  that  the  frequency  response 
between  response  at  i’  and  excitation  at  k  is  the  same  as  the 
function  between  response  at  k  and  excitation  at  i.  Equation  19 


is  frequently  written  in  the  form: 

2?T  „  K 


r  .r* 

ih _  +  ik 


where 


r=l  (s  -  sr)  r=l  (s-sr)  (s-s*) 

AT,  =  Residue  at  pole  s  (i.e.^i  ^  k  ) 
xx  r  a 


The  impulse  response  of  the  system  can  be  obtained  from  Equation 
21  by  performing  an  inverse  transform  to  obtain: 


**ik  ' "}  ~ 


Kik  e‘ 
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I 


Since  the  roots  occur  in  conjugate  pairs i  Equation  22  can  be 


written  in  the  form 


Hi*")’  2 


4|  e'  Sr  “r1  cos[(«r^l-  J/MIk]  (23) 


where 

4  -  <*L 

Equation  2'j  indicates  that  the  impulse  response  of  the  system 
can  be  represented  by  a  summation  of  the  number  of  damped 
cosine  waves  times  the  appropriate  modal  parameters. 

The  multi -degree-cf-freedom  (i.'iDCF)  curve  fitting  procedures 
in  the  nodal  analysis  program  calculates  the  value  of 
in  the  above  equations.  Therefore,  in  the  case  where  aF^ 
was  determined  from  a  displacenent/f orce  frequency  response 
function,  the  value  of  ar  can  be  determined  from  the  equation! 

SpF  vj£ 


where 

aF,  is  determined  from  a  displacement/force  function. 

X  X 


In  the  case  where  a  velocity/force  frequency  response  function 

was  curve  fit  with  the  MDOF  orccedure,  the  rarameter  a_  is 

*  r 

determined  from  the  following* 


*1 _ 


x  j  CO. 


where*  cJ„  is  in  trie  units  of  rad/sec  ar.c 

AFk  is  determined  from  a  velocity/force  function 
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Similarly,  if  an  acceleration/force  frequency  response  function 
is  used,  ar  i?  determined  from  the  equation* 

,  2> 


ar  = 


x  (-«/) 


where 

<Or  is  in  units  of  rad/sec  and 

Af,  is  determined  from  an  acceleration/force  function 

If  an  analytical  model  is  to  be  created  from  the  test  data, 
the  parameters  a^,  'i'r,  £Jr  and  are  all  that  is  necessary 

to  describe  the  component  with  complex  normal  modes.  However, 
in  some  cases  an  analyst  would  like  to  use  a  "real"  mode 
approximation  with  the  associated  effective  mass  or  effective 
stiffness  in  order  to  describe  the  component  under  test  via 
the  following  equation* 

^  vl/  *•"  vl/r 

H.,  =  _ Z* _  (24) 

4  lk  rarE  H2*  W2  ~  ^rw"r] 

In  that  case  it  is  recommended  that  this  approximate  representa¬ 
tion  be  determined  by  setting  the  magnitude  of  the  mode  shape 
coefficient  equal  to  the  magnitude  of  the  complex  mode  shape, 
value  and  the  sign  of  the  mode  shape  coefficient  from  one  of 
the  following  procedures* 

1)  Inverse  of  the  sign  of  the  imaginary  portion  of  the  mod 
shape  coei.fici.ent  when  a  displacement/force  freouency 
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response  function  is  used  to  determine  the  node  shape 
coefficients. 

2)  The  sign  of  the  real  portion  of  the  node  shape  coeffi¬ 
cient  when  a  velocity/force  frequency  response  function 
is  used  to  define  the  node  shape  coefficients. 

3)  Sign  of  the  imaginary  portion  of  the  node  shape 
coefficient  when  an  acceleration/force  frequency 
response  function  is  used  to  determine  the  node  shape 
coefficients. 


The  effective  mass  necessary  to  approximate  the  actual  frequency 
response  with  that  described  by  Equation  24  can  be  determined 
from,  one  of  the  following  equations: 


— 

r 


(Approx.  'J’f)  (Approx. 
where  A?.  is  determined  from  a  dis?lacenent/f orce  func 


nr  = 


(Approx,  )  (Approx.  \s  v) 

’*  K*l  ' 

where  is  determined  from  a  veiocity/force  function 


mr  " 


(Aoorcx.  ^i)(  Aourox.^k) 


Alk 


x  cj 


M  ^ 

where  is  determined  from  a  acceleratior./forcc  fur.c 


Ir.  order  to  represent  a  component  in  an  overall  system  model 
via  a  "real”  node  approximation,  the  following  approach  can 
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..frequently  be  used#  The  uncoupled  equations  of  motion  for 
the  component  in  terms  of  modal  coordinates  [  %  ]  aret 

£-to2['m/}  +  jco["cO  +  [*]  *  C?i  3 

where  nu]  is  a  diagonal  matrix  of  effective  masses 

['  c*]  is  a  diagonal  matrix  of  effective  dampness 

[%k„]  is  a  diagonal  matrix  of  effective  stiffnesses 

The  motion  of  the  physical  coordinates  [q]  is  related  to  the 
motion  of  the  modal  coordinates  byi 

[Q]  «['*']  [O 

Symbolically,  this  can  be  represented  by  the  following  diagram 


Structural 

Coordinates 

Modal 

Coordinates 


*1 


Therefore,  a  component  can  be  represented  analytically  from  test  data 
in  an  overall  system  model  by  a  set* of  springs,  masses,  dampers  and 
equations  of  constraint  which  relate  the  motion  of  the  physical  coor¬ 
dinates  to  the  motion  of  the  modal  coordinates.  Since  the  equations 
of  constraint  can  be  quite  volumnous,  the  KASTaA.M  input  and  :--ATSIX 
Seneration  task  provides  the  capability  to  generate  KASrSA;*  kuiti  ?oi 
Constraint  (I.FC)  equations  in  a  relatively  automatic  manner. 
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The  frequency  response  in  a  specified  range  can  be  approx- 
iaately  described  in  terns  of  the  following  quantities i 

1)  “Residual  Inertar.ee"  of  the  nodes  of  vibration 

below  the  range  of  interest. 

2)  The  nodes  of  vibration  which  are  resonant  in  the 
specified  frequency  range. 

3)  "Residual  Compliance"  of  the  modes  of 'vibration 
above  the  range  of  interest. 

Mathematically  this  can  be  expressed  ast 


“lie  *  #  + 


lill 


r=l  ar(Jw+^rur  -  <3wr  A- 


This  concept  is  shown  graphically  in  the  following  figure si 
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The  frequency  response  in  a  specified  range  can  be  approx¬ 
imately  described  in  terms  of  the  following  quantities » 

1)  "Residual  Inertar.ee"  of  the  modes  of  vibration 

below  the  range  of  interest. 

2)  The  modes  of  vibration  which  are  resonant  in  the 
specified  frequency  range. 

3)  "Residual  Compliance"  of  the  modes  of  vibration 
above  the  range  of  interest. 

Mathematically  this  can  be  expressed  asi 

••  =  its  .  #•  ¥£  _ .  , 

"in  u1  T  jSr  ar(j **$rur  *  j<jr  A-sJ‘)  T  'ik 

This  concept  is  shown  graphically  in  the  following  figures i 


MSwjaiKiniain 


j  Iw-inertar.c; 

e 


FrW«YjM,)  |  Ff#cu#nev  ^  | 

of  Intiftn 

The  3DC?  ar.d  M3C?  curve  fitting  procedures  available  in  the 
Estimation  task  are  used  to  evaluate  the  contribution  due  rc 
the  modes  ’which  are  resonant  in  the  frequency  range  under 
investigation.  In  order  to  determine  the  ccncri’cucicn  of  the 
residual  effects  the  Generate  Residual  command  in  the  Frequency 
Response  Synthesis  task  is  used. 
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APPENDIX  C 

Mechanical  Analysis  Data 


THOUGH 


S-G3: 


appendix  d 

Basic  Program  plus  Results 


Program  Variables 


D2  :  Logarithmic  decrement 


f  s  Natural  frequency 


S  :  Section  Modulus  (ln^) 


-  Ll. l, 

VTT? 


v  -  2  rad/sec  for  CORT 

-  3  rad/sec  for  RYERSON 

■  2.5  rad/sec  for  BLOUGH 


SM  -  96169  for  CORT 

-  42196  for  RYERSON 

-  60000  for  BLOUGH 


D3  s  Damped  Period  ■  tj  ■  6 _ 

wn  *  * 


°P-T 

-  BM  -  SM  *  2 

12  *  2000 


M  :  Bending  Moment 


-D3- 

10  DIM  DU00,4),T*64,D*<100)64 
30  INPUT  L*,fJ 
30  FOR  1=1  TO  N 

40  READ  D*<I),0(I,2>,D(1,3),D(1,4) 

SO  NEXT  1 

60  PRINTUblNC  61 

61  X  LOCATOR  CODE  DAMPING  FREQ  LOG.  DEC  DAMPED  PERIOD  B.M. 

62  PRINTUSING  63 

63  2  FACTOR  (HZ)  (SEC)  (FT -TONS) 

70  FOR  L*1  TO  N 

80  M=2 

100  IF  D(L,t1)=0  THEN  155 

110  DE^ttPl^DiL.fD/d-tDtL.rDtH)^^) 

120  P=E 

121  3=92169 

122  IF  L<19  1HHN  130 

123  P=3 

124  S =42196 

1E5  IF  L<38  THEN  130 

126  p =2 . 6 

127  3=60000 

130  D3=D2/(F*U(L,M) ) 

136  M*  (S*  (DO. ,4  ) /2 ) )  /  (£000*12 ) 

140  PRINT  USING  141,G¥(L)  ,D<L,2)  ,D(L,3>  ,02  ,D3,i1 

141  7.  *##0006664  #.#####6  6.6#  6.66666  666.66#  6.666tm 

ISO  COTO  160 

156  PRINTUSING  1S6,D*<L> 

166  7.  6«#66t5##66  .  .  .  . 


160  NEXT  L 
170  END 

1B0  DATA  ‘CHSL* 18*20",. 018,. 29, 14666 

190  DATA  "CH5L -‘21*33*  ,  .0126  ,.33,0.,  "CHSL  *20*23’-  , .  0246  ,  .43 , 12236  ,  "0:53*15*02*  , .  01 

25..  34. 2 0306 

200  DA) A  "CH3B*1S*22* ,.0131, .34, 16464, "CHS£*1S»24*" , . 0145, .3, 20036, rCM5L*EC»03" , 

.  00804,.  22, 9160,  “CMS'_*20*01" ,  0.0, 0.0,0. ,  "C MEL* 19* OS" .  .0CE2S..2E  ,55  £4 

210  DATA  "CM3»<1S*3J"  ,  .0167  ,  .35 , 10395 , *0)56* 15*06  *  ,0.0, 0.0,0. , ‘CHIB* 17*27 " , .005) 

1 . . 32 .9317 , *CLSL* 17*31* , .0137 , .22 ,3925 , *CL5L*0C*1£* , . 0173 , .SB ,3925 , "C uSc* 02*21 ’ , 

.00719,. £2,4335  ' 

220  LATA  "CLSE«09*10‘  ,  .0122,  .31 ,5179 , *CLSB* 10*07 ‘ ,  .0203 ,  .31 ,43S£ ,  "CL5E-*  19*25"  ,.0 

151. . 3.4595 

221  DATA  *  KHi-*5-l*  1  ’,.00423,.  57,1 1800,  "RHL*4  -2*1  *,.00455,.  56, 9 110,  THE  *5-1*12“  , . 
00433 , .57 ,13420 

222  DATA  "RHfc»b-l*13*  ,  .00514 , .57 , 12621 , *RHB*5-1* 15*  ,  .00573  ,  .57 ,13960 , "R-  '■•'»5-l*i6 
",.00522,. 57, 11800, *WIL*5-2* 19" , .00539, .56, 7550, *RML*5-£ '24* ,.00745, .55, /730, ' RM 
L*5-2*39" , .00395, .57 ,7090 

223  DATA  *KM2*4-1*10“,. 00*05,. 53, 6000, *RME-4-l*ll" ,.00362,. 58, 5900, "RMB*4 -1*15" , 
.00337, .57,6540," RnB*  4  - 1* 17  * , . 0  0475 , . 57 . 794  0 

£24  DATA  "RLL* 6-1*3' , .00625 , .57 ,4470 , "RLL*6- 1* 9* , . 00495 , .56 ,4390 ,  *RLL*6-1*10  ' , . 0 
0524,. 44, 4300 

225  DATA  "RLB* 1-1*21" ,  .00359  ,  .58  ,4170, "RLB* 1-1*24* ,  .00335 , .59 ,3730 , ■RLB#1-1*00" , 
0,0,0 

22to  DATA  "BHL* 13*51" ,0,0,0, "DHL* 13*50* ,. 0171 , .323 , 10330 , *E<HB* v8R*36" , . 0174 , .42 ,2 
0026 , "BUB* 138*45" , . 0232 , .39 , 162  74 , * 5H3*  0ER*3E  * , . 033 , .39 , 14225 

227  .DATA  "BML*13*S2* ,  .0303 ,  .382 ,7581 ,  *  BALMS' 53"  .  .0242  ,  .379, 5971, "EML* 13*47'  ,.01 
2B,. 337, 7050  ,*tMB«13*3*,.  0143,. 41, 10721,  *E;’.3 ‘13*5*,.  0327,. 4, 10C03,*B:.£*  13*4  “,.01 

2. . 41.10111 

222  DATA  "BLL*  13 '46*  ,  .0202 ,  .38 ,4  V:  9  , "BLL*13*43*  ,0,0 ,0 ,  TLL»i3«45" ,  .0105  .  .363  ,27'- 
7 ,'  BlB*SR*28*  ,  0,0,v),*BLD»8R*3)J"  , .  0117,.  43, 4042,  "BLB>i3P.- 57"  ,.0164,.  4^,4  211 
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FACTOR 

(HZ) 

s 

(SEC) 

(FT -TONS] 

CHSL* 13*20 

0.018000 

0.29 

0.11517 

3.199 

2.816E+04 

CHSL *21*33 

0.012500 

0.33 

0.07953 

3.181 

O.OOOE+OO 

CHSL *2 0*28 

0 . 024600 

0.43 

0.15846 

3.220 

2.349E+04 

CHSB*15*02 

0.012500 

0.34 

0.07553 

3.181 

3.899E+04 

CHSB*15*££ 

0.018100 

0.34 

0 .11582 

3.199 

3.161E+04 

CHSB*19*£4 

0.014600 

0.30 

0.09309 

3.188 

3.C47E+04 

CHSL«20*03 

CMSL*£0*01 

0.003640 

0.22 

0.05092 

3.167 

1.756E*04 

0.008280 

0.28 

0.05245 

3.167 

1.834E+04 

CMSL*19*0S 

CMSB*15*3D 

C(1SB*1S*06 

0.016700 

0.35 

0.10671 

3.194 

2.092E+04 

0.006110 

0.32 

0.03S62 

3^160 

1.88SE+04 

CMSB*1?*£7 

CLSL*17*31 

0.013700 

0.28 

0.08727 

3.185 

7.651E+03 

CLSL*06*12 

0.017300 

0.22 

0.11061 

3.196 

7 .555E+03 

CLSL*08V'21 

0.007190 

0.22 

0.04550 

3.164 

8 .324E+03 

CLSB*09*10 

0.012200 

0.31 

0.07760 

3.180 

9.944E+03 

CLS3*10*07 

0.020300 

oisi 

0.13019 

3.206 

8.414E+03 

CLSB*19v-2b 

0.015100 

0.30 

0.05633 

3.189 

B.823E+03 

RHL*b-l*l 

C .004830 

0.5/ 

0.03045 

2.104 

1.037E+04 

RHL«4-2*1 

0 .004560 

.0.56 

0.02871 

2.103 

8 .003E+03 

RH6*5-1*12 

0.004330 

0.57 

0.02732 

2.103 

1 . 179E  +  C4 

RHB*S-1* 13 

0.00S140 

0.57 

0.03246 

2.105 

1.109E+04 

RHB*5-1*15 

0.005730 

0.57 

0.03552 

2.106 

1 .227E  >04 

RHB*S-l*lfe 

0.005220 

0.57 

0.03297 

2.105 

1 . 037E<  04 

RML*b-2*19 

0.005390 

0 .56 

0.03404 

2.105 

6.637E+03 

KnL*b-2*c4 

0.007450 

0.56 

0.04716 

2.110 

6 .7955  +  03 

RML*S-£*39 

0.003950 

0.57 

0.02491 

2.102 

6.232E+03 

RKB»4-1*10 

0.004050 

0.58 

0.02555 

2.102 

5 .274E  +  03 

RM6*4-l*li 

0.003620 

0.58 

0.02282 

2.102 

5 . 1S6H+03 

Hrib*4-l*15 

0.003870 

0.57 

0.02441 

2.102 

5.749E+03 

Rr.B*4-l*17 

0.004750 

0.57 

0.02393 

2.104 

6.979E+03 

RLL*5-1*8 

0.006250 

0.57 

0.03951 

2.107 

3.925E+03 

RLL*6-1*9 

0.004950 

0.56 

0.03125 

2.104 

3.859E+03 

RLL»6-1*10 

0.005240 

0.44 

0.03309 

2.105 

3 .780E  +  03 

RLB* 1-1*21 

0.003590 

0.58 

0.02253 

2.101 

3.665E+03 

HLB* 1-1*24 
RLB* 1-1* 00 

0.003360 

0.59 

0.02118 

2.101 

3.322E+03 

BHL*13*bl 

0.017100 

0.32 

0.10931 

2.556 

1.353E+04 

BHU* 13*50 

BHB*0b‘R*3S 

0.017400 

0.42 

0.11126 

2.557 

2.S03E+04 

8HB*13R*4b 

0.023200 

0.39 

0.14923 

2.572 

2 . 034E*  04 

BHB* 032*35 

0 . 033000 

0.39 

0.21442 

2.599 

1.778E+04 

BML* 13*52 

0.030300 

0.38 

0 . 19632 

2.591 

9.476E*03 

UML* 13*53 

0.024800 

0.37 

0.15978 

2.577 

7 .463E  +  C3 

BML*13*47 

0.012800 

0.32 

0.08146 

2.545 

8 .862E  +  03 

BfiB*13*3 

0.014300 

0.41 

0.09115 

2.549 

1 .340E+04 

BH8*13*5 

0.022700 

0.40 

0.14594 

2.571 

1 .325E  +  04 

6PIB*13*4 

0.012000 

0.41 

0.07631 

2.543 

1 .263E+C4 

BLL* 13*45 

0.020200 

0.38 

0.12553 

2.565 

5.573E+03 

BLL* 13*43 

0. 010300 

0.38 

0.06539 

2.539 

4. 72 IE *03 

BLL* 13*45 

BLB*8R*£8 

0.011700 

0.43 

0.07433 

2.543 

5.052E+C3 

BLB*8R*38 

BU>*13K'S? 

0.018400 

0.42 

0.11777 

2.560 

5 . 138E  +  03 

j.niia 


appendix  e 

American  Bureau  of  Shipping  letter 


JLr  <3M  &  roooff 

.'/tyrU  YNC/bv 

/&/ 

23  January  1979 


Lt.  Mark  Noll 

Department  of  Transportation 
U.S. Coast  Guard  (G-DSA-1/TP44 ) 

Washington,  D.C.  20590 

Dear  Mark: 

We  have  received  the  package  containing  a  sample  copy  for 
the  purpose  of  transfer  of  fund  and  your  report  "Evaluation 
of  SDRC  Damping  Analysis"  dated  18  October  1978.  Matters 
regarding  the  feasibility  study  are  now  being  handled  by 
Dr.  D.  Liu,  our  Chief  Engineer.  I  believe  he  has  discussed 
this  matter  with  you  on  the  phone  already. 

I  would  like  to  offer  some  comments,  as  you  so  requested,  on 
the  SDRC  work  and  your  report  as  follows.  I  believe  that 
the  SDRC  work  and  our  efforts  {hereby  abbrieviation  as  ABS ' 
work  for  easy  reference)  on  analyzing  the  damping  coefficient 
through  the  use  of  sea  trial  data  (stress  history)  are  based 
on  the  same  principle,  a  happy  coincidence.  There  are,  however, 
major  differences  as  noted  below: 

(1)  SDRC  uses  a  multi-degree-of-freedom  system  and  ABS  uses  a  one 
d.o.f.  system.  It  is  believed  that  the  inter-modal  cor¬ 
relation  in  the  power  spectrum  is  neglibible  with  the  provision 
of  (2)  as  a  safeguard. 

(2)  We  employ  partial  spectra  based  on  a  concept  proposed 
by  Vanmarcke  (OTC  1971)  rather  than  the  full  spectrum. 

(3)  ABS'  work  also  incorporate  the  equilibrium  branch  of 

a  wave  spectrum  which  is  assumed  to  be  inversely  proportional 
to  the  fifth  power  of  frequency.  This  is  also  a  departure 
from  Vanmarcke' s  original  idea. 

(4)  Since  the  power  spectrum  transformed  from  the  stress 
time-history  is  expressed  as  blocks  with  constant  ordinate 
within  a  frequency  cell  (bin) ,  the  fabricated  spectrum  should 
also  be  in  the  same  form  as  in  the  case  in  our  analysis.  Our 
matching  criterion,  as  proposed  by  Vanmarcke,  is  to  have  the 
quantity  p  of  the  fabricated  spectrum  to  match  that  of  the 
measured  spectrum  where 

p  -  1  -  m, 2/ (m  m  ) 

1  0  2 
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TO  Lt.  Mark.  Noll  DA,E  1/23/79  mE  *tF 

Department  of  Transportation 
U.S.  Coast  Guard 


and 


m^  =  ith  spectral  moment. 

In  pursuing  the  curve  fitting,  the  damping  coefficient  as  well 
as  the  resonance  frequency  (not  necessarily  equal  to  the 
central  value  of  the  frequency  cell)  are  regarded  as  variables. 
In  doing  so,  the  two  spectra  will  have  the  same  zero,  first 
and  second  spectral  moments. 

We  have  employed  542  intervals  of  the  "Cort"  data  set  (1973 
season)  from  which  one  value  of  the  damping  coefficient  is 
obtained  for  each  interval  and  the  resulting  histogram  is 
shown  in  Fig.  1. 

It  should  be  noted  that  not  all  results  are  reliable.  Some 
intervals  have  spectral  ordinates  so  small  that  they  are 
nothing  but  small  ripples  in  the  power  spectrum  (no  springing) . 
For  this  reason,  cutoff  values  for  the  RMS  value  (square  root 
of  m0)  have  been  used  and  the  resulting  mean  and  standard 
deviation  of  the  set  of  results  are  shown  as  functions  of  this 
cutoff  value  as  displayed  in  Fig.  2.  As  a  result,  those 
unreliable  intervals  are  screened  out  as  the  cutoff  RMS  value 
increases  and  both  the  mean  and  standard  deviation  exhibit 
stablized  trend.  In  this  figure,  the  curves  labeled  as 
"C  £  0.039"  represent  further  screening  of  elimininating 
the  unusually  high  values  of  the  damping  coefficients  which 
amounts  to  less  than  5%  of  all  intervals  computed.  It  is 
evident  that  the  damping  coefficient  is  virtually  independent 
of  the  stress  level  as  the  "stablized"  portion  of  the  curves 
indicate.  It  suffices  to  apply  this  conclusion  without  limit¬ 
ing  it  to  the  high  stress  samples  since  only  the  cases  of  high 
springing  stress  are  of  interest. 

It  is  apparent  although  the  number  of  intervals  considered  in 
the  SDRC  calculation  is  very  small  compared  with  our  sample 
size,  that  the  mean  value  of  their  results  as  well  as  the 
spread  agree  very  well  with  our  results  (i.e.,  mean  *  0.0145, 
standard  deviation  =  0.0045). 

As  indicated  in  your  report,  results  obtained  utilizing  the 
logarithmic  decrement  method  are  much  higher  than  those  from 
spectral  synthesis.  We  have  observed  the  same  tendency  through 
the  case  of  the  RANDOMDEC  method  which  yields  a  vibration 
signature  in  the  time  domain.  Such  time  domain  signatures 
are  subsequently  analyzed  by  way  of  logarithmic  decrement. 

One  possible  explanation  of  the  discrepancy  between  results 
obtained  in  these  two  methods  is  that  in  the  time  domain,  the 
signal  is  not  truly  monochromatic.  In  this  regard,  perhaps 
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Lt.  Mark  Noll 

Department  of  Transportation 
U.S.  Coast  Guard 


DATE 1/23/79  FILE  REF 


the  test  of  a  ship  with  zero  surge  would  yield  closer  agree¬ 
ment  between  the  two  approaches . 

Based  on  the  discussion  presented  above,  it  appears  that  the 
spectral  synthesis  approach  is  more  promising.  In  this  regard, 
the  approximate  shape  of  the  wave  energy  spectrum  should  also 
be  accounted  for  as  our  analysis  indicates  that  the  omission 
of  this  aspect  leads  to  about  20%  higher  values  for  the  damping 
coefficient  on  the  average. 

I  hope  the  preceding  comments  would  provoke  more  subsequent 
discussions  on  this  subject.  If  you  have  any  suggestions 
or  comments,  please  feel  free  to  contact  me. 


Enel . 


Sincerely  yours, 


Chen 
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APPENDIX  F 

Power  Spectral  Densities  submitted 
Jay  SDRC 
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